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Recent experiments show that spin-polarized current may influence the state of generally accessory 
element of spin- valves - an antiferromagnetic (AFM) layer, which is used for "pinning". Here we 
study the dynamics of AFM component of the "pinned" ferromagnetic (FM) layer induced by 
simultaneous application of the spin-polarized current and external magnetic field. 

We find stability range of such a configuration of FM/AFM system in which orientation of FM 
magnetization is parallel to AFM vector. We calculate the field dependence of the critical current 
for different orientations of the external magnetic field with respect to the crystal axes of FM/AFM 
bilayer. We show the possibility of stable current-induced precession of AFM vector around FM 
magnetization with the frequency that linearly depends on the bias current. Furthermore, we 
estimate an optimal duration of the current pulse required for switching between different states of 
FM/AFM system and calculate the current and field dependencies of switching time. The results 
obtained reveal the difference between dynamics of ferro- and antiferromagnets subjected to spin 
transfer torques. 

PACS numbers: 72.25.-b, 72.25.Mk, 75.50.Ee 
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I. INTRODUCTION 

Spin transfer torque (STT) is the torque that is applied 
by non-equilibrium spin-polarizcd conduction electrons 
onto a magnetic layer'" '. This effect creates, for exam- 
ple, the ability to switch nanoscale magnetic devices at 
GHz frequencies and stimulate emission of microwaves by 
steady electric current'-^'. Key elements of the spintronic 
devices, that enable an information coding, control and 
manipulation by an electric current, are two ferromag- 
netic (FM) layers. The "pinned" layer acts as a polarizer 
for conduction electrons, while the state of a "free" layer 
may be altered by STT. However, recent experiments''"^" 
give an indirect evidence that the spin-polarized current 
may also influence the state of another, generally acces- 
sory, element - an antiferromagnetic (AFM) layer, which 
is used for "pinning" . 

The characteristic value of current density at which 
the switching of AFM state takes place varies from 
10^ A/cm^ (for an insulating AFM and current-in-plane 
geometry'') to 10^ A/cm^ (for a metallic AFM and 
current-perpendicular-to-plane geometry'') and, in prin- 
ciple, can be smaller than the critical current density 
in the similar giant magnetoresistive structures without 
an AFM layer (oc 10^ ^ 10^ A/cm^, Refs. 11-13). Spin- 
polarized current also affects both the exchange bias and 
the coercive field of "free" FM layer ' ' . Combined ap- 
plication of spin-polarized current and external magnetic 
field gives rise to various switching scenarios depending 
on the thickness, sequence and material of FM and AFM 
layers^ ^. On the other hand, the physical mechanism and 
details of such a nontrivial dynamics are still unclear. 



Due to the efforts of several theoretical groups^''"^*^ the 
concept of STT is extended to the systems with i) dif- 
ferent types of magnetic ordering, including nonuniform 
and disordered FM (that, in principle, could be further 
extended to the magnetic systems with noncoUinear and, 
probably, AFM spin ordering); ii) different nature of the 
magnetic ordering and interaction between the charge 
carriers and spins, i.e., sd-exchange in the magnets with 
localized spins or itinerant magnetism in the transition 
metals. It also became clear that the STT phenomena 
could result from the atomic scale spin dependent scat- 
tering (i.e., hopping of a conduction electron between the 
sites with different directions of the magnetic moments). 
In some particular cases an AFM can even work as a 
polarizer for conduction electrons and exert spin torque 
on the adjacent FM or AFM layer, as it was predicted 
in Rcfs.lii and 17. However, all of the published calcu- 
lations are based on assumption of quantum coherence 
and so, are applicable to the perfect samples with ideal 
interfaces. 

In our previous paper^" we proposed the phcnomcno- 
logical model that describes the current-induced phenom- 
ena in AFMs on the same footing as in FM materials. It 
was assumed that the total angular momentum is con- 
served during an interaction of spin-polarized transport 
electrons with each of magnetic sublattices. 

In the present paper we apply this model for the de- 
scription of the precessional switching processes induced 
by simultaneous application of the spin-polarized current 
and external magnetic field to an AFM component of the 
"pinned " layer depicted in Fig. 1(a). Our chief aim is to 
study the different static and dynamic regimes of AFM 
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layer and to find the way to induce a stable precession of 
an AFM vector starting from a certain configuration of 
FM/AFM bilayer. We also try to find similar and differ- 
ent features in the current-induced dynamics of FM /FM 
and FM/AFM bilayers. We anticipate our approach to 
be a starting point for a more comprehensive analysis 
of the multilayered magnetic systems in the presence of 
high-density current. For example, joint behavior of the 
FM and AFM layers could be analyzed, with account of 
the exchange bias coupling. 



II. SPIN TRANSFER TORQUE IN THE 
MULTISUBLATTICE MAGNETS 

According to Berger^ and Slonczewski^' the physi- 
cal mechanism of STT in fcrromagnets can be explained 
in the following way. When a free electron transverses 
(or reflects from) an interface between the nonmagnetic 
(NM) and FM layers, its spin state can be reversed due 
to exchange interaction with the localized magnetic mo- 
ments of FM. This process results in rotation of the lo- 
calized moments in a way that ensures conservation of 
the total spin of the system. 

Generalization of the Berger's and Slonczewski's ideas 
to the case of multisublattice materials is not straightfor- 
ward due to more complex character of the magnetic or- 
dering. Particularly, in AFMs the direction of the atomic 
magnetic moments varies on the length scale of atomic 
distances leading to zero net magnetization if averaged 
over few lattice constants. However, just as in FMs, the 
spin-polarized electrons transfer spin torques on each of 
atomic sites"''!^'-"'^!. 

The magnetic structure of AFM may be described with 
the use of a few macroscopic vectors Mj (in the simplest 
case j = 1,2) called the sublattice magnetizations (per 
unit volume) that are formed due to strong exchange 
coupling. So, it seems reasonable to assume that while 
entering an AFM, the conduction electron transfers spin 
angular momentum to any of the magnetic sublattices 
(see Fig. lb). Corresponding STT Tj exerted by the 
j-th sublattice is then presented in a standard form as 
follows: 



M, 



[Mj X [Mj X pcur]], 



(1) 



where J is the current spin-polarized in Pcur direction, 
|Pcur| = 1, the constant aj = eh'y /{2MojVe) is propor- 
tional to the efficiency e of scattering processes, V is the 
volume of AFM region, h is the Plank constant, e is the 
electron charge, 7 is the modulus of the gyromagnetic 
ratio, and Mg^ = |Mj| is the saturation magnetization of 
j-th sublattice (the value of Moj is supposed to be un- 
changed under external fields). Positive current (J > 0) 
corresponds to injection of electrons into AFM layer. 

Then, the dynamics of AFM can be described by a 
set of Landau-Lifshitz- Gilbert equations for Mj vectors 



supplemented with the Slonczewski term (1): 



-7[M, X H,] 



M, 



[M, X M, 



"^^■"^■[Mj- X [Mj X peur]], 
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(2) 



where H , 



-dw/dM.j is the "generalized force" (an 
effective local field acting on the magnetic moment of a 
sublattice) and w is free energy (per unit volume) of an 
AFM layer. For the sake of clarity we describe relaxation 
of an AFM layer in the simplest form, with the use of a 
single Gilbert damping parameter aQ equal for all mag- 
netic sublattices (although the relaxation mechanisms in 
AFM crystals are very complicated and diverse"^.) 

The last two terms in r.h.s. of Eq. (2) are responsible 
for dissipation processes in the AFM layer. To illustrate 
this fact we calculate the rate of free energy losses in 
assumption that dissipation is small and in zero approx- 
imation Mj = -7[Mj X Hj]. Thus, 



(3) 



E 



'^3 J 



Mf - — ( Pcur • [Mj X Mj 



In principle, Eqs. (2), (3) could be used for descrip- 
tion of different complicated magnetic structures (com- 
pensated AFMs, weak FMs, ferrimagnets) . In the limit- 
ing case of the completely equivalent sublattices (Mi ~ 
M2 = • . .) the set of equations (2) turns into a standard 
Landau-Lifshitz-Gilbert-Slonczewski equation for FMs. 

In the particular case of AFM with two magnetic sub- 
lattices it is more suitable to rewrite Eqs. (2) in terms of 
macroscopic magnetization (FM vector) m = Mi + M2 
and AFM order parameter (AFM vector) 1 = Mi - M2: 



m = 7([HA/ xm] + [HL x 1]) 



2Mn 



m X m 



[ixi] 



crJ 



([m X [m X ]] + [1 X [1 X Pcur]]) 



2Mo 

i = 7([HM X 1] + [Hl X m]) 

crJ 



(4) 



2Mo 



([m X [1 X pc 



2Mo 

[1 X [m X Pcur]]) 



[m X 1] + [1 X m] 



(5) 



Here Hm = —dw/dm is an effective magnetic field 
within an AFM layer that includes an external magnetic 
field. Hi — ^dw/dl is a magnetic anisotropy field conju- 
gated to an AFM order parameter, and |Mi| — IM2I = 
Mo. 

Equations (4), (5) describe the dynamics of FM and 
AFM vectors in the presence of spin-polarized current 
and generalize the Landau-Lifshitz-Gilbert-Slonczewski 
equation for the systems with more than one magnetic 
sublattice. In what follows we base our considerations on 
these equations. 
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III. MODEL 

Let us consider a pinned layer (Fig. la) of a typi- 
cal exchange-bias spin- valve that includes an AFM layer 
whose thickness cJafm is much smaller than the char- 
acteristic scale of the magnetic inhomogeneity. On the 
other hand, ^afm is large enough to ensure an AFM 
ordering within the layer. High current densities are 
achieved in a small region (10-^100 nm in diameter), in 
which both FM and AFM layers could be considered as 
a single domain. In the case of a moderate pinning (i.e., 
when the magnetic anisotropy of FM is comparable with 
the unidirectional anisotropy induced by exchange bias) 
the FM works as a spin-polarizer whose state is not af- 
fected by the precession of AFM vector in the adjacent 
layer. So, magnetization of the FM layer is assumed to 
be fixed and is described by the vector Pcui - 




Figure 1. (Color online) Effect of spin transfer torque 
within a pinned layer, (a) General structure of a pinned 
layer. Due to exchange coupling, an AFM layer pins the direc- 
tion of magnetization pcur of the adjacent FM layer. In turns, 
the FM layer polarizes the spin current flowing to AFM. (b) 
Transfer of spin torques Ti,2 from free electrons (e~) to sub- 
lattice magnetizations Mi and M2 (solid arrows - before and 
dotted arrows - after an interaction with the conduction elec- 
trons). STT produces small FM moment m which, in turn, 
sets up a solid-like rotation of Mi, M2 around m (arc ar- 
rows), according to Eq. (6). (c), (d) Sketch of eigen modes 
of an AFM vector for different orientations of the external 
magnetic field. Easy axes (EA) are parallel to the film plane. 

In our analysis we take into account the fact that 
AFMs (e.g., FcMn, IrMn) widely used in spintronic de- 
vices show strong exchange coupling (corresponding ex- 
change field He ^ H^) between the different magnetic 
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sublattices that keeps magnetizations Mi and M2 al- 
most antiparallel even in the presence of the external field 
Ho <C He- In this case the state of AFM is described by 
the only vector order parameter 1, and far below the Neel 
temperature |1| « 2Mq. Spin-polarized current and/or 
external magnetic field induce small tilt of the sublattice 
magnetizations (Fig. lb) formally described by the FM 
vector m. Vector |m| ^ |1| is a slave variable and can 
be expressed from Eq. (5) as follows (see Appendix A for 
details): 



1 



[i X 1] 

2-fMoHE 2MoHe 



[1 X [Ho X 1]], 



(6) 



where Hq is the external magnetic field (which, particu- 
larly, can be induced by Oersted field of a current). 

Substitution of the expression (6) into Eq. (4) gives 
rise to a closed equation for AFM vector: 



[ixl]=7{(2i(Ho,l)-[lx[Hoxl]] 
-7[Ho x l](Ho, 1) + 27Moif£[Hi x 1] 

+aGHE[\ X i] + aJHE[\ X [1 X Pcur]]} 



(7) 



Equation (7) describes a solid- like motion of AFM vector 
in which |1| is almost unchanged. Nevertheless, due to ad- 
ditional (compared to FM) degrees of freedom, this equa- 
tion differs from the standard Landau-Lifshitz-Gilbert- 
Slonczcwski equation for FM vector M (see (2) for 
in order (includes the 2-nd order derivative of the dy- 
namic variable, 1 instead of M). 

Analysis of the last term in (7) shows that the current- 
induced contribution is proportional to the AFM order 
parameter (vector 1). Thus, one can expect the influence 
of spin-polarized current on the dynamics of AFM vec- 
tor to be at least as strong as in FM materials (other 
things being equal). Moreover, the current-dependent 
(STT) term in the r.h.s. of Eq. (7) contains large mul- 
tiplier He- This is manifestation of the so-called effect 
of exchange enhancement when some interactions (e.g., 
gap values, spin-phonon coupling) are more pronounced 
in AFMs than the analogous interactions in FMs'^'^. 



IV. DYNAMICS OF AFM VECTOR WITHIN 
THE LAGRANGE APPROACH 

An effective formalism for investigation of AFM dy- 
namics is based on the use of Lagrange formalism^ 
Equation (7) can be regarded as an Euler-Lagrange equa- 
tion of the second kind in the presence of dissipative ex- 
ternal forces (see Appendix B). Corresponding Lagrange 
function has a form 



1 



FM 



-1^ 



1 



Aj^MoHe 2jMoHe 



Ho • [1 X i] 



1 



4.M0HE 



[1 X Ho] 



.(I)- 



(8) 
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Here WanO) is the energy of magnetic anisotropy (per 
unit volume). 

To take into account the effect of STT that can work 
both as a source or drain of energy for an AFM layer, we 
deduce from (3) the dissipative Rayleigh function (see 
Appendix B) 

HafM = T^i' - fpcur ■ [1 X []) (9) 



(10) 



47M0 2jMo 
that describes the rate of the energy losses 



dw 



1 



d\ 



Analysis of dissipative function (9) shows that STT 
phenomena in AFM have one general property which is 
not peculiar to FM. While STT always changes the en- 
ergy of FM layer, some types of motions in AFM could 
be nondissipative even in the presence of spin-polarized 
current. Linearly polarized oscillations of the vector 1, 
sketched in Fig. Ic, give an example of nondissipative 
mode (neglecting the internal damping) . And, vice versa, 
the most effective energy pumping induced by the cur- 
rent takes place for any precessional, circular polarized 
motion of AFM vector in the plane perpendicular to the 
direction of current polarization Pcur (see Fig. Id). 



V. STABILITY DIAGRAM 

To illustrate the peculiarities of non-dissipative and 
dissipative current-induced dynamics, we analyze sta- 
bility of the state with parallel orientation of AFM 
and FM vectors, Pcur||l, for two different configurations 
of the external magnetic field Hq, depicted schemati- 
cally in Figs. lc,d. For the definiteness. an AFM layer 
is supposed to have slightly tetragonal (almost cubic) 
anisotropy induced, e.g., by the shape effects or/and in- 
teraction with the neighboring layers (including possible 
influence of the exchange bias). Two easy axes {Z and 
Y) are parallel to the film plane. In this case the mag- 
netic anisotropy energy is modeled with the following 
expression 



Mn 



X 



8A/3 



(4+4+4), (11) 



where i/an|| is the intrinsic anisotropy field within the film 
plane and the small out-of-plane anisotropy field iJan_L ^ 
^an|| is responsible for weak tetragonality of the sample. 

In the absence of field and current a single AFM 
layer has two equivalent equilibrium orientations of AFM 
vector (see Fig. 2a, b): \\\Z and (as can be eas- 
ily obtained from minimization of the magnetic energy 
(11)). Correspondingly, a FM/AFM bilayer has two 
stable configurations- ' with l||pcur (Fig. 2a) and with 
1 _L Pcur (Fig. 2b) . These two configurations should have 
different macroscopic properties (e.g., different magne- 
toresistance, different exchange bias field, etc.) and in 



this sense are analogous to the parallel (P) and antipar- 
allel (AP) configurations of FM/FM multilayers (Fig. 2 
c,d). In analogy with FM/FM systems, the reversible 
switching between the l||pcur and 1 + Pcur states can be 
achieved by application of the external magnetic field to 
the free (in our case, AFM) layer. 

The switching field should be oriented parallel to AFM 
vector, its critical value coincides with the spin-fiop tran- 
sition field Hs-i = 2^ H,^^\\He for AFM layer (also ex- 
change enhanced). 



A Z 





Figure 2. (Color online) Stable magnetic configurations 
witli different macroscopic properties. In FM/AFM bi- 
layer (a,b) the magnetization (||pcur) of FM layer is fixed, 
AFM vector can be switched from (a) parallel to (b) perpen- 
dicular orientation with respect to pcur by application of the 
external magnetic field Hq > Hs-f within an AFM layer. In 
FM/FM bilayer (c,d) the magnetization of free layer can be 
switched from (c) parallel to (d) antiparallel orientation with 
respect to pcur by the external magnetic field Ho applied an- 
tiparallel to Pcur- In both cases (c) and (d) the switching can 
be also induced by current. 

When the current is injected into bilayer, configuration 
with Pcur II 1 is still equilibrium, but not necessarily stable. 
To find the values of critical current and field, we analyze 
the frequencies of eigen modes of AFM layer (magnetiza- 
tion of FM layer Pcur and, correspondingly, the current 
polarization being fixed). 



A. Configuration Ho -L 1 

In the crossed initial orientation Hq + 1 (Fig.lc) the 
linearized equations of motion for the generalized coordi- 
nates Ix, W {Iz « 2Mo - (4 + 4)/2) take the following 
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form: 



Ix + 27AFM^X 



j%) Ix + iHeoJW = 0, 



ly + 27AFM^y + ^^yW - "/HectJIx = 0, 



(12) 



Here 7afm = iHeolg/'^' ^ ^x,y is a damping co- 
efficient tliat can be estimated from the line-width of 
AFM resonance, ujh = 7^0- The values ujx = 
2jy/{Han± + H^^\\)He, ujy = 2'-fy/H^^i\HE are the eigen 
frequencies of free oscillations in the absence of field and 
current that also could be measured in AFM resonance 
experiments. It is worth noting that the values of eigen 
frequencies are enhanced due to exchange coupling (mul- 
tiple He) compared to analogous values in FM with the 
same value of anisotropy field. 

Equations (12) describe the case when the magnetic 
field is directed along the hard anisotropy axis, Ho||X. 
Configuration with HoUy (field is parallel to an easy 
axis) is treated in an analogous way. 

It can be easily seen from (12) that below the critical 
current 



1 



2jHe(7 



I 2 2 , 2 I 

\uJx - t^y + UJij\. 



(13) 



the eigen modes have linear polarization and corre- 
spond to oscillations of vector 1 within XZ or YZ plane 
(Fig. Ic). In this case the spin torque transferred from 
the current affects the eigen frequencies of spin excita- 
tions, 

1 



^±=2 +(^Y+^H 



I I 2 2 

± lUJx — ^ 




(1) 



(14) 



but does not affect the effective damping coefficients (as 
it is the case in FM). 

It should be stressed that in the absence of external 
field the value of critical current depends upon anisotropy 



{uj^ — bjy) oc i?an± of thc magnetic interactions within 
and perpendicular to the film plane. The magnetic field 
applied perpendicular to the vector 1 enhances (if Hq is 
parallel to an easy axis) or weakens (if Hq is parallel to 
a hard magnetic axis of AFM) the effective anisotropy. 
So, magnetic field can be used for control of the critical 
current. If anisotropy is weak (77anJ_ ^ ^an||j or, cquiv- 
alcntly, \ujx — wy| ^ ^y), it can be effectively reduced 
with the field whose value is much less than thc spin-flop 
one. Ho < ifs-f • 

Above the critical current, |J| > Jcr\ polarization of 
free oscillations changes from linear to circular (elliptic) 
and STT contributes into the energy dissipation. For 
one of two modes of free oscillations the current-induced 
pumping competes with thc internal damping. Starting 
from the critical value 



4? ^ J(j. 



272 



AFM 



{CO 



'HI 



(15) 



the average energy losses per oscillation period are neg- 
ative (pumping is greater than damping) and the state 
with the parallel alignment of current polarization and 
AFM vector becomes unstable. 

f 2) 

As seen from Eq.(15), the value of critical current Jd- 
is independent on the directions of current, field and spin 
polarization (Pcurr), in contrast to the threshold current 
for FM. 



B. Configuration Ho||l 

Another type of dynamics is observed in the case when 
the field is applied parallel to 1 (Fig. Id). In this case 
polarization of eigen modes is circular (or elliptic) even 
for J = 0, as follows from symmetry considerations and 
from analysis of equations of motion written in terms of 
appropriate generalized coordinates l± = Ix ^ iW'- 



J 



l± + 2 (7AFM T i^^n) l± 



(16) 



So, oscillations of 1 can actively take up an energy from 
the current and STT affects the damping coefficient, not 
the frequency of oscillations. Instability point is attained 
as soon as the spin-polarized current overcomes the effect 
of internal friction. 

Fig. 3a) shows the field-current stability diagram for 
the case of isotropic AFM (iJani = and, correspond- 
ingly, ujx = ^y)- Within the shaded area 



J 



27AFM^fo 



aHp 



< 



27AFM-ffs-f 



\Hn\<H.. 



(17) 



the static state with l||pcur| 
value, |J| > |Jcr|, where 



Jr. 



27. 



AFM 



Z is stable. Above the critical 



crHf 



{Ho + Hs-isignJ) , 



(18) 



the current may keep up a stable rotation of AFM vector 
around Pcur (Fig- 3b). Sign reversal of STT (resulted 
from the reversal of either direction of current, J ^ — J, 
or direction of polarization, Pcur —> — Pcur) gives rise to 
rotation in opposite direction. If the field value is greater 
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Figure 3. (Color online) Dynamics of AFM vector in- 
duced by steady spin-polarized current, (a) Stability 
diagram under combined action of field and current. Static 
states of the bilayer (schematically shown with an arrow for 
Pcur and double arrow for AFM vector) are stable within the 
shaded area. Strong current (|J| > |Jcr|) induces precession 
of an AFM vector around the current polarization Pcur, di- 
rection of rotation depends on the sign of J. Strong field 
{\Ho\ > Hs-i) induces spin-flop transition combined with the 
current-induced precession, (b) 3-D evolution of AFM vector 
1 in overcritical regime ( J > Ja) in the presence of magnetic 
field Hq. In the initial state the AFM vector is slightly mis- 
aligned from Z axis. Under the action of STT the vector 1 
spirals away from Z axis with a steadily-increasing preces- 
sion angle and the angular frequency Q. Final state corre- 
sponds to precession with the stable frequency within XY 
plane {Iz = 0). (c), (d) Time dependence of the angular fre- 
quency fl and Iz projection. Arrows indicate the moment, 
tmin, at which monotonic decrease of Iz switches to the de- 
caying oscillations. Envelope (dash line) corresponds to re- 
laxation (oc exp(— 7APMi)) caused by internal damping. 

than spin-flop field. \Hq\ > i?s-f, the state with l||Ho 
is unstable even in the absence of current and in the 
final state the AFM vector is perpendicular to Hq and 
Pcur- These results also keep true for small but nonzero 
anisotropy i?an±- 

The current-induced precession of 1 vector is also sta- 
ble in the "high-field" region, |_ffo| > ^^s-f- However, 
detailed analysis of the dynamical phases and transition 
lines in this region is out of scope of this paper. 

Some features of the current-induced instability in con- 
figuration with l||Ho arc similar to those observed in 
FM/FM bilayers. First, in both cases the stability region 
is defined by an internal friction which stands up against 
the current-induced rotations-''. Second, the value of 
critical current linearly depends on the field^"'-'. Thus, 
application of the magnetic field results in variation of 
the critical current and opens a possibility to reduce Jcr, 
as seen from Fig. 3a. 

On the other hand, there are still few principal dif- 
ferences between FM/FM and FM/AFM bilayers listed 
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(a) 




Figure 4. (Color online) Switching between the differ- 
ent configurations of FM/FM (a) and FM/AFM (b) 
bilayers. Magnetization of the fixed layer is shown with ma- 
genta (thick) arrow, that of the free layer with violet (thin 
arrow), double arrow shows orientation of AFM vector, (a) 
Switching between P and AP states can be achieved by the 
field or current applied in two opposite directions, (b) Tran- 
sition from parallel to perpendicular configuration can be in- 
duced by current (arbitrary direction) and field applied along 
initial orientation of AFM vector. Transition from perpendic- 
ular to parallel configuration can be induced by the field only. 



below. 

• In the FM/FM bilayer (one FM layer being fixed) 
switching between P and AP states can be achieved 
by application of either field or current (Fig. 4a). 
In contrast, in the FM/AFM bilayer switching be- 
tween l||pcur and 1 ± Pcur states can be achieved 
by a combined application of field and current 
(Fig. 4b). Namely, the current induces transition 
only from l||pcur to 1 _L Pcur because the last state is 
stable in the presence of current. To reverse 1 vec- 
tor back to l||Pcur configuration, one needs to apply 
an external field Hq > Hs-i parallel to 1 (spin-flop 
transition). 

• In the FM/FM bilayer the direction of current 
(from fixed to free layer or opposite) is important, 
P — > AP and AP P transitions take place at op- 
posite directions of current (Fig, 4a). In contrast, 
in the FM/AFM bilayer destabilization of l||pcur 
state takes place irrespective of the current direc- 
tion (Fig. 4b). However, an external magnetic field 
removes such a degeneracy. 

• The bilayers with AFM should show exchange re- 
duction of the critical current compared to FM/FM 
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bilayers providing that the free FM and AFM lay- 
ers have the same magnetic resonance frequencies 
(or anisotropy field of FM is close to spin-fiop field 
of AFM) and the same quality factor (= cj/7afm), 
as can be seen from Eq. (18). 

VI. DYNAMICS IN OVERCRITICAL REGIME 

A FM layer subjected to the direct spin-polarized cur- 
rent shows one interesting effect - stable precession of 



magnetization with the angular frequency close to the 
frequency of spin- wave mode'' '. To find out whether 
such an effect could be observed in AFM, we consider in 
details the dynamics of AFM vector in overcritical regime 
(I J| > I Jcrl) assuming that Pcur||Ho||Z. 

We use the standard parametrization of AFM vector 
with the spherical angles 6 and Lp, Ix = 2A/o sin0cos(p, 
ly = 2A/osin(?sini^, Iz = 2MqCOsO, to deduce the fol- 
lowing dynamic equations: 



+ 27AFM^ + sin cos 

— \{lp - loh) sin^ d] + sin^ 
at 



{ip - loh) + 



2 , '^2^JtV(^^^Qg2^)- '^^ '^g sin^g(7-f-cos4^) 



0, 



sin 2lp 



— s\v? 9 sin 



(19) 



As it was already mentioned, an AFM under considera- 
tion is an oscillator with the high quality factor {lox.y ^ 
7afm)- In other words, energy dissipation takes place on 
the time scale much greater than the characteristic pe- 
riod of free oscillations. In this case for analytical treat- 
ment of Eqs. (19) one can apply the asymptotic method 
of rapidly rotating phase originated by Bogolyubov and 
Mitropolskii^''* . 

According to this method, the motion of AFM vector is 
decomposed into rapid rotation with the frequency 17 cx 
u!x,Y and slow variation of amplitude and frequency with 
the characteristic time scale oc 1/7afm- In the simplest 
case of isotropic AFM (i?an± = 0, or ux ~ loy) the 
only rapid variable is <^ = ^[t)t. Equations for slow 
variables Q.{t) and 9{t) (J7, 9^0) are obtained from 
(19) by averaging over the period of rotation: 



+ 27afm6' 
-t- sin 9 cos 9 
d 



•J 

cJq — (r2 — ujh) — —loq sin^ ' 



dt 



= 0, (20) 

[{n ~ ujh) sin^ 9] + sin^ 9 - iHeJcf) = 0. 



If, in addition, Cl ^ 9, the first of equations (20) describes 
ID motion (dynamic variable 9) in a potential well (see 
Fig. 5) 



u{9-n) 



1 



1 



2 ■ 2 

LOf^ sm ( 



(21) 

with the friction defined by coefficient 7afm • The second 
of Eqs. (20) describes the current-induced variation of 
both variables 9 and fi. 

Equations (20) have two interesting solutions. The first 
one, corresponds to the circular polarized free oscillations 
of AFM vector with an amplitude 6* = 6'o <C 1 and eigen 
frequency f^o = oJx +i^h- However, in overcritical regime 



an amplitude 9o growth with an increment proportional 
to the offset from the critical current value: 



7AFM 



Jc: 



(22) 



The second solution with 9 = Tr/2 corresponds to 
steady rotation of AFM vector in XY plane {Iz ~ 0) with 
the angular frequency Sloo = {J/ Jcy)^o- Energy dissipa- 
tion per period of rotation is zero, due to the pretty bal- 
ance between the magnetic damping and current-induced 
pumping. This solution is stable when \J\ > | Jcr|, as can 
be seen from analysis of the potential U{9; fl). Small de- 
viations of AFM vectors from XY plane (16* - 7r/2| < 1) 
relax due to internal friction as 



cos (7 oc e 



— 7AFMi 



cos(f2ei -I- ip), 



(23) 



where ip is a, phase that depends upon initial conditions 
and 



ne = \lni^-2n^Lu„ + -{ojl 



'Hi 



(24) 



So, the state of steady precession is approached during 
the time 1/7afm that depends upon the internal mag- 
netic damping. 

To illustrate all the described peculiarities of AFM dy- 
namics in the presence of spin-polarized current we solve 
the original Eqs. (19) numerically with the initial condi- 
tions 9 = 9q = 0.001, LP = 7r/2, = and = f^o. In 
other words, at i = an AFM vector deflects from equi- 
librium orientation \\\Z through the small angle 9q <^ 1 
within the ZY plane. Initial velocity corresponds to 
that mode of free oscillations which is unstable for the 
chosen current direction. For calculations wc used the 
following dimensionless values: lox = 'j-'y = 6.28 and 
7AFM = 0.314 (that corresponds to the quality factor 
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Figure 5. Profile of the effective potential well (21), arbitrary 
units. l-f2<f2o,2-n>no. Arrows indicate static solution 
(Q — 0) unstable in overcritical regime and stable stationary 
solution {9 = 7r/2). 

20). Time unit equals to the period of free oscillations in 
the absence of field and current. 

Fig. 3b shows a typical trajectory of AFM vector (nor- 
malized to a unit length) in the presence of steady current 
J = 2.5Jcr < and field H = 0.2Hs-{. With the de- 
scribed initial conditions, the motion of 1 vector starts as 
a rotation around Z axis with the eigen frequency ^Iq of 
free oscillations. Due to the energy pumping from STT, 
an amplitude of oscillations (1 projection on XY plane) 
slowly increases with an increment r (see Eq. (22)). 

The final state {t oo) corresponds to the above 
described steady rotation of AFM vector in XY plane 
(Iz = 0) with the angular frequency floo- In analogy 
with FM, such a precessional state of an AFM layer can 
be a source of spin waves. In contrast to FM, the an- 
gular frequency fioo is proportional to the current value. 
The absolute value of fioo is greater than the charac- 
teristic spin- wave frequency (oc J7o) which in AFM can 
range THz values (e.g., for bulk FeMn the energy gap is 
7 meV"'' that corresponds to linear frequency cx2 THz). 
So, FM/AFM bilayer can be considered as a potential 
emitter of high frequency spin waves. 

Figs. 3c and 3d illustrate the time evolution of the rota- 
tion frequency Q and the component Iz between the ini- 
tial and final states. Due to nonlinear effects, defiection 
of 1 from the initial direction is accompanied by decrease 
of 17. At a certain moment t = imin (shown by arrow 
in Figs. 3 c,d) a monotonic decrease of Iz changes into 
decaying oscillations around the average value Iz ~ 0- 
Relaxation of 1 to XY plane is due to internal damping 
and follows the law cx exp(— 7AFMi) (see Eq. (23) and 
envelope in Fig. 3d) . The details of relaxation to the pre- 
cessional state are shown in Fig. 6 where we compare an 
exact solution (points) of (19) with the asymptotic form 




' 1 ' 1 ' 1 

25 30 35 40 

t 



Figure 6. (Color online) Relaxation dynamics of slow ampli- 
tude Iz a.t t > tmin- Points - numerical simulation, solid line 
- approximation according to Eq. (23), dashed line - envelope 
oc exp(-7AFMi). 

(solid line) calculated from the expressions (22), (23). 

In our simulations of AFM dynamics we found that not 
only Izj but the rotation frequency il, energy dissipation 
rate and the effective potential energy averaged over a 
period of rotation (= 27r/17) has an extrcmum at the 
moment t ~ imin- This means that the system passes 
through the crossing between two attraction points in the 
phase space. So, we interpret the time interval tmin as a 
switching time between two stable states of the FM/AFM 
bilayer (see Figs. 2 a,b, and 4 b). The exact value of 



r 100 




2 I -1 ■ 1 ' l-fr 
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cr s-f 

Figure 7. (Color online) Switching time from parallel 
to perpendicular configurations shown in Fig. 2 a,b. 

(a) Current dependence of the calculated switching time, tmin, 
and the effective pumping coefficient r oc \ J — Jci\~^ , external 
magnetic field Ho = 0.2Hs{ is parallel to Z. Inset shows corre- 
lation between r and tmin. (b) Field dependence of switching 
time tmin at fixed bias current (points). Solid line shows ap- 
proximation oc |_H"sf -I- -ffo|~^. 
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imin depends upon the initial deflection Oq of AFM vector 
from Z axis (or, in other words, from the amplitude of 
spontaneous fluctuations and, hence, from temperature). 
However, the current and field behaviour of imin is the 
same for different initial conditions and correlates with 
the current (Fig. 7a) and field (Fig. 7b) behavior of the 
characteristic time r of dcstabilization (22). Though the 
absolute values of tmin and r are different, they show 
good correlation (see inset in Fig. 7a) in rather wide range 
of current values. Both values decrease moving far from 
the stability point This opens a way to diminish 

the switching time by increasing the current value or by 
decreasing the critical current with the magnetic field. 



,/ 


a 







■ 10 20 30 40 50 





10 20 30 40 50 



VII. CONTROLABLE SWITCHING OF AFM 
STATE 

From practical point of view it is important to achieve 
a controllable switching between the different equilibrium 
states of the FM/AFM bilayer (say, l||pcur ^ 1 -L Pcur or 
1|| ~ Pcur) using the current pulses of minimal duration 
and amplitude. We investigate dynamics of AFM vec- 
tor under the rectangular current pulses (schematically 
shown by red line in Fig. 8) of different duration and am- 
plitude in overcritical regime (| J| > | Jcr|)- If pulse dura- 
tion is below iniin, AFM vector returns back to its initial 
state after the current is switched off. 

Fig. 8 demonstrates the switching processes initiated 
by the current pulse J = 2.5 Jd- with the duration slightly 
greater than imin- The chosen pulse duration ensures 
maximum deflection of AFM vector from the initial di- 
rection as seen from time dependence of Iz (Fig- 8 a). 
After the current is switched off, the AFM vector re- 
laxes to the final static state within XY plane through 
the damping oscillations during the time 1/7afm- Reg- 
ular rotation around Z axis supported by current also 
vanishes with the end of current pulse, as seen from f2 
behaviour (Fig. 8b). The final orientation of AFM vec- 
tor is parallel to one of the easy axes within XY plane 
(in the presence of field, Fig. 8c) or can be also anitipar- 
allel to the initial 1 direction (180° switching to Z easy 
axis) in the absence of external field (Fig. 8d). Due to 
degeneracy, the final state is very sensitive to the initial 
conditions and pulse duration and can be predicted only 
statistically. 



VIII. CONCLUSIONS 

In summary, we studied the dynamics of AFM layer in 
the presence of spin-polarized current and external mag- 
netic field. 

On the basis of a simple model of slightly tetrago- 
nal AFM with two magnetic sublattices we demonstrated 
the following features of current-induced behavior in the 
FM/AFM bilayer. 



Figure 8. (Color online) Dynamics of AFM vector in- 
duced by a pulse of spin-polarized current, (a) Time 
dependence of Zz (thick or blue line) under the rectangular 
current pulse (thin or red line). Pulse duration close to tmin 
is enough to ensure maximal deflection of Iz from the initial 
value, (b) High frequency rotation of AFM vector persists as 
long as current induces STT, as seen from time dependence 
of il. Just after the current is switched off the frequency of 
rotation goes down to zero and AFM vector lies down to XY 
plane. Relaxation time (= 1/7apm) of both SI and Iz is due 
to internal damping, (c) Magnetic field applied parallel to the 
initial orientation of 1 induces 90° switching, the final state of 
AFM vector is parallel to XY plane, (d) In isotropic AFM in 
the absence of field a spin-polarized current may induce 180° 
switching. 



1. Spin transfer torque induces the loss of stability in 
the FM/AFM system if AFM vector is parallel to 
FM magnetization. Configuration with AFM vec- 
tor perpendicular to FM magnetization is stable 
in the presence of current and field (in accordance 
with the prediction made in Ref.17). This means 
that the high density current can induce reorienta- 
tion of AFM vector if an angle between 1 and Pcur 
differs from 7r/2. 

2. Such a current can also induce a stable precession 
of AFM vector in the plane perpendicular to FM 
magnetization. The frequency of precession is of 
the order of frequency of free oscillations and lin- 
early depends on current. 

3. The value of critical current can be tuned by ap- 
plication of the external magnetic field. The value 
of switching time can be tuned by both field and 
current. 

We anticipate the same features in noncoUinear AFM 
with three (like IrMua and Mn3NiN) and more (like 
FeMn) magnetic sublattices. 

On the other hand, our model predicts irreversible 
current-induced switching between the parallel and per- 
pendicular orientations of 1 and Pcur- This result 
is the consequence of the chosen tetragonal magnetic 
anisotropy. In the case when an angle between easy 
axes of AFM differs from 7r/2 (hke in NiO or FcMn), 
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the current-induced switching seems to be possible be- 
tween all the configurations of 1 and Pcur as long as 1 has 
a nonzero projection on Pcur- 

The described response of AFM vector to electron cur- 
rent may change the properties of the pinned layer. Due 
to the weak but nonzero exchange coupling between AFM 
and FM layers, reorientation or precession of AFM vec- 
tor results in variation of the exchange bias field and, 
consequently, gives rise to the shift of switching fields of 
spin- valve. We illustrate this effect qualitatively in Fig. 9. 
Consider a typical spin-valve with the pinned FM layer. 
Suppose, an AFM layer is inhomogeneous (multidomain) 
and high density current gives rise to reorientation of 
AFM vector in some of domains (Fig. 9a). The ratio of 
the rotated domains is proportional to the integral cur- 
rent. On the other hand, reorientation of some AFM 
domains results in diminishing of the exchange bias field 
that keeps FM magnetization of the pinned layer. Varia- 
tion of the bias field is also proportional to the ratio of the 
rotated AFM domains. So, in the presence of current the 
critical field at which magnetization of the pinned layer 
is reversed decreases linearly. However, linear shift of the 
bias field can be also induced by STT between FM layers. 
If AFM layer is not affected by spin transfer torque, the 
stability region of AP configuration increases for one cur- 
rent direction and diminishes for an opposite, as shown 
in Fig. 9b. On the contrary, if spin torque is transferred 
to an AFM layer and is not transferred between two FM 
layers, the stability region of antiparallel configuration of 
FM layers diminishes for any current direction (Fig. 9c). 
The last type of current dependence (among others) was 
observed in the experiments Ref.l4. Decrease of the ex- 
change bias field irrespective of current direction was also 
observed in Rcf.!). 

Linear shift of the bias field induced by the current was 
observed in nanopillars ' that included coupled permalloy 
(FM) and FeMn (AFM) layers. In these experiments 
a combined application of the magnetic field and high- 
density current resulted in an increase of the exchange 
bias field from -100 to 100 Oc. 

Another evidence of STT effects in AFM could be 
found from the detailed analysis of the field/current 
dependence of magnetoresistance, as it was done in 
Rcfs.l), 8-10. Magnetoresistance of spin valve should de- 
pend on the angle between FM and AFM vectors (in ad- 
dition to the dependence from mutual orientation of FM 
vectors in "free" and "pinned" layers) and can change 
due to the current-induced switching of AFM vector. 
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Figure 9. (Color online) Current dependence of the ex- 
change bias field, (a) Exchange bias spin-valve structure 
with multidomain AFM layer. AFM vector in some of do- 
mains deflects from the initial orientation due to STT. (b) 
The shift of bias for typical spin transfer (AFM layer is not 
afi'ected by current), (c) The shift of bias in the case when 
spin torque is transferred to AFM layer. 



where He is the spin-flip field of the exchange nature, 
Wan is anisotropy energy, Hq is external magnetic field. 
In assumption that HeMq ^ (or, equivalcntly. 

He 3> Hl) and He ^ Hq (strong exchange couphng), 
the effective magnetic field is expressed as follows: 



H 



M 



dw 
dm 



Hn 



He 

2Mn 



(A. 2) 



Wc substitute (A. 2) into 
(with ac and J) and get 



neglect dissipative terms 



1 = 7 [Ho X 1] 



He 
2A/o 



[m X 1] 



(A. 3) 



To obtain an explicit expression (6) for m we multiply 
both sides of Eq.(A. 3) by the vector 1 and take into 
account that [1 x [m x 1] = mF « AMqUI. The last re- 
lation follows from the fact that below the Neel temper- 
ature both vectors 1 and m are bound by the constraint 
P + m' 

requirement |Mi| = IM2I = Mq 



= 4A/q, (1, m) = (that is equivalent to the 



Appendix A: Relation between magnetization and 
AFM vector 

Small macroscopic magnetization |m| <C |1| of AFM 
layer can be excluded from Eq.(5) in the following 
way'^"^''*". Free energy of AFM layer is modeled as 

He 

W = TTT™^ + ^^an " Hom, (A. 1) 



Appendix B: Lagrange and Rayleigh functions for 
ant ifer romagnet 

The Lagrange function (8) and the Rayleigh function 
(9) are selected in such a way that to fulfill the following 
requirements: 

i) The dynamic Eqs.(7) are the Euler-Lagrange equa- 
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tions of the second kind with dissipative forces: 



losses (3) according to Eq.(lO). 



FM 



X 1 



d OCa: 
dt di. 

di. 



BCa 



FM 



X 1 



X 1 



(A. 1) 



ii) The effective potential energy in Lagrange function 
coincides with the magnetic anisotropy energy (11). 

in) Rayleigh function is related with the rate of energy 



Expression for energy losses is obtained from (3) by sub- 
stitution Ml = — M2 = 1/2. Contribution from m- 
depending terms into Rayleigh function is small and so, 
is neglected. 

Lagrange approach makes it possible to account for 
the constraint |1| = 2Mq (valid far below the Neel point) 
by appropriate choice of two generalized coordinates qk 
{k ~ 1,2) instead of three components of vector 1, as 
described in the paper. 



1 L. Berger, Phys. Rev. B 54, 9353 (1996) 

^ J. C. Slonczewski, Phys. Rev. B 39, 6995 (1989) . 

^ J. Slonczewski, J. Mag. Mag. Mater. 159, LI (1996). 

* S. Kaka, M. R. Pufall, W. H. Rippard, T. J. Silva, S. E. 
Russek, and J. A. Katine, Nature 437, 389 (2005). 

^ A. Slavin and V. Tiberkevich, IEEE Trans. Magn. 45, 1875 
(2009) 

® Z. Wei, A. Sharma, A. S. Nunez, P. M. Haney, R. A. Duine, 
J. Bass, A. H. MacDonald, and M. Tsoi, Phys. Rev. Lett. 
98, 116603 (2007). 

^ S. Urazhdin and N. Anthony, Phys. Rev. Lett. 99, 046602 
(2007) 

* X.-L. Tang, H.-W. Zhang, H. Su, Z.-Y. Zhong, and Y.-L. 
Jing, Appl. Phys. Lett. 91, 122504 (2007) 

^ N. V. Dai, N. C. Thuan, L. V. Hong, N. X. Phuc, Y. P. Lee, 
S. A. Wolf, and D. N. H. Nam, Phys. Rev. B 77, 132406 
(2008) . 

^° J. Bass, A. Sharma, Z. Wei, and M. Tsoi, Jour, of Mag- 
netics (Korean Magnetics Society) 13, 1 (2008) 

" M. Tsoi, A. G. M. Jansen, J. Bass, W.-C. Chiang, M. Seek, 
V. Tsoi, and P. Wyder, Phys. Rev. Lett. 80, 4281 (1998). 
J. Grollier, V. Cros, H. Jaffres, A. Hamzic, J. M. George, 
G. Faini, J. Ben Youssef, H. Le Gall, and A. Pert, Pliys. 
Rev. B 67, 174402 (2003) 

D. C. Ralph and M. D. Stiles, J. Mag. Mag. Mater. 320, 
1190 (2008) 

J. Basset, A. Sharma, Z. Wei, J. Bass, and M. Tsoi, Proc. 

of SPIE, Spintronics 7036, 703605 (2008) . 

Y. Tserkovnyak, H. J. Skadsem, A. Brataas, and G. E. W. 

Bauer, Phys. Rev. B 74, 144405 (2006) 
^® A. Nunez, R. Duine, P. Haney, and A. MacDonald, Phys. 

Rev. B 73, 214426 (2006). 
" P. M. Haney and A. H. MacDonald, Phys. Rev. Lett. 100, 

196801 (2008) . 

1* Y. Xu, S. Wang, and K. Xia, Phys. Rev. Lett. 100, 226602 
(2008) 



E. V. Gomonay and V. M. Loktev, Sov. J. Low Temp. 
Phys. 34, 198 (2008) . 
^° P. M. Haney, D. Waldron, R. A. Duine, A. S. N. nez, 
H. Guo, and A. H. MacDonald, Phys. Rev. B 75, 174428 
(2007) . 

P. M. Haney, D. Waldron, R. A. Duine, A. S. N. nez, 
H. Guo, and A. H. MacDonald, Phys. Rev. B 76, 024404 
(2007) 

^2 V. G. Bar'yakhtar, Sov. Phys. — JETP 60, 863 (1984) 
A. L Akhiezer, V. G. Bar'yakhtar, and S. V. Peletminskii, 
Spin Waves, Interscience (Wiley) ed., North- Holland Se- 
ries in Low Temperature Physics, Vol. 1 (North-Holland, 
Amsterdam, 1968). 

L Bar'yakhtar and B. Ivanov, Sov. — J. Low Temp. Phys. 
5, 361 (1979). 

In the bilayer structure the orientation of AFM vector is 
governed by competition between the intrinsic magnetic 
anisotropy (11) and exchange coupling with the adjacent 
FM layer. However, in the typical cases an effective field of 
exchange coupling is much smaller than the spin-flop field 
of AFM and thus can be neglected. 

S. I. Kiselev, J. C. Sankey, I. N. Krivorotov, N. C. Emley, 
R. J. Schoelkopf, R. A. Buhrman, and D. C. Ralph, Nature 
425, 380 (2003). 

Y. V. Gulyaev, P. E. Zil'berman, E. M. Epshtein, and R. J. 
Elliott, Sov. Phys.— JETP 100, 1005 (2005). 
N. N. Bogolyubov and Y. A. Mitropolskii, Asymptotic 
methods in the theory of nonlinear oscillations (Taylor and 
Francis, Inc., 1961) 548 p. 

Y. Endoh, G. Shirane, Y. Ishikawa, and K. Tajima, Solid 
State Commun. 13, 1179 (1973). 
^° A. V. Kimel, B. A. Ivanov, R. V. Pisarev, P. A. Usachev, 
A. Kirilyuk, and T. Rasing, Nature Physics 5, 727 (2009) 



